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Chapter  1 
INTRODUCTION 


This  technical  note  describes  a  simplified,  computer  efficient.  Wind  Factor 
Simulation  Model  (WFSM)  developed  by  the  United  States  Air  Force  Environmental 
Technical  Applications  Center  (USAFETAC)  in  response  to  a  request  for  such  a  model  bv 
the  Military  Airlift  Command  (MAC).  MAC  desired  this  model  for  inclusion  in  their 
large  scale  airlift  system  simulation  effort,  COLOSSUS. 

The  WFSM  calculates  climatological  wind  factors  by  Sawyer's  equivalent  headwind 
techniaue  (Air  Weather  Service,  1977)  for  arbitrary  great  circle  routes  at  specified 
altitudes  for  any  of  four  seasons  of  the  year  in  any  of  three  modes  or  wind  options 
(calm  wind,  mean  wind,  and  90-oercent  worst  wind).  In  its  present  form  the  WFSM  does 
not  provide  climb  or  descent  winds.  Nor  does  it  provide  temperatures  and  aircraft 
performance/fuel  consumption  information.  The  simulated  wind  option,  intended  to 
provide  variability  in  the  winds,  has  not  been  implemented,  and  requests  for  the 
simulated  wind  are  treated  as  requests  for  the  mean  winds. 

The  Wind  Factor  Simulation  Model  consists  of  seven  FORTRAN  subprograms  totaling 
approximately  2500  lines  of  source  code,  including  extensive  comments.  The  package 
includes  subprograms  for  calculating  great  circle  distance  and  initial  great  circle 
headina.  Procedure  plus  data  for  the  seven  subprograms  occupy  5. IK  words  of  core 
storage  on  a  36-bit  general  purpose  computer.  The  WFSM  assumes  that  a  wind  data  base 
has  been  read  from  a  tape  or  disk  file  into  a  COMMON  block.  Hence,  core  storage  for 
a  file  control  block  and  buffer  must  be  added  to  the  above  total  to  obtain  a 
realistic  estimate  of  the  core  storage  required  to  run  the  wind  factor  model. 
Obtaining  a  wind  factor  for  a  long  route  requires  approximtely  0.1-0. 5  seconds 
central  processing  unit  (CPU)  time  on  a  Honeywell  6080  general  purpose  computer. 
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Chapter  2 

WIND  FACTOR  SIMULATION  REQUIREMENTS 


2 . 1  General 


Requirements  stipulated  for  the  Wind  Factor  Simulation  Model  (WFSM)  were 
originally  quite  extensive.  Negotiations  with  the  user  reduced  this  list  of 
requirements  to  those  stated  below.  In  the  future,  a  number  of  the  oriqinal 
requirements  removed  by  negotiation  are  expected  to  reappear.  Anticipated  future 
requirements  are  discussed  in  section  2.5  below. 

2.2  Required  Outputs 

The  wfsm  must  produce  simulated  route-mean  wind  factors  for  reasonably  arbitrary 
great  circle  routes  between  point  "A"  (takeoff)  and  Doint  "B"  (landing)  at  either  of 
two  constant  altitudes  (25,000  ft  and  35,000  ft)  for  any  of  the  four  seasons  of  the 
year  (winter,  spring,  summer,  and  fall)  in  any  of  three  modes  or  wind  options:  (a) 
calm  wind,  (b)  mean  wind,  and  (c)  90-percent  worst  wind.  A  fourth  mode,  the 
simulated  wind  case,  is  one  of  the  future  requirements  discussed  in  section  2.5 
below.  If  simulated  winds  are  requested  by  the  user,  mean  winds  will  be  provided. 

Specifically,  the  Wind  Factor  Simulation  Model  must  output  the  route-mean  qround 
speed  in  knots  (kt)  for  a  given  route-mean  airspeed  in  knots  and  calculated  route- 
mean  wind  factor  in  knots.  In  addition,  the  model  must  calculate  the  qreat  circle 
distance  (GCD)  in  nautical  miles  (NM)  between  any  two  points  whose  global 
(latitude/longitude)  coordinates  are  known. 

A  capability  possessed  by  the  model  because  of  internal  needs,  but  not 
specifically  requested  by  the  user,  is  the  ability  to  calculate  the  initial  heading 
in  degrees  along  a  great  circle  route  between  any  points  whose  global  coordinates  are 
known.  Table  1  provides  a  general  summary  of  WFSM  outputs.  More  specific  output 
descriptions  are  provided  in  Chapter  5. 


TABLE  1.  SUMMARY  OF  MODEL  OUTPUTS. 


Reouired/Internally  Needed  Quantity  Units 

Route-mean  ground  speed  kt 

Great  circle  distance  (GCD)  NM 

Initial  great  circle  heading  dec 


2.3  Inputs  Provided 


A  user  employing  the  WFSM  to  calculate  a  route-mean  wind  factor  and  ground  speed 
must  stipulate  the  global  (latitude/longitude)  coordinates  of  points  "A”  (takeoff) 
and  "B"  (landing),  the  Julian  base  date  (1-366)  of  the  wind  factor  reauest  (used  to 
determine  season  of  the  year),  aircraft  altitude  (25,000  ft  or  35,000  ft),  mode  or 
wind  option  (calm,  mean,  or  90_percent  worst)  and  airspeed  in  knots. 

To  accommodate  future  growth  the  user  must  also  stipulate  Greenwich  mean  base 
time  (0.0-24.0)  of  the  wind  factor  request  and  forecast  hours  ahead  of  the  wind 
factor  required.*  Presently,  these  inputs  are  ignored  and  may  have  any  value. 
Table  2  provides  a  general  summary  of  WFSM  inputs.  More  specific  descriptions  are 
provided  in  Chanter  5. 


•Greenwich  mean  base  time  is  the  Greenwich  mean  time  for  which  the  wind  factor  is 
desired.  In  the  present  WFSM,  Greenwich  mean  base  time  is  the  same  as  Greenwich  mean 
time.  Future,  more  sophicated  versions  of  the  WFSM  will  produce  wind  factor  fore¬ 
casts.  Then  the  Greenwich  mean  base  time  will  be  the  Greenwich  mean  base  time  for 
which  the  wind  factor  forecast  is  desired  minus  the  number  of  forecast  hours  ahead. 
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TABLE  2 


SUMMARY  OP  MODEL  INPUTS. 


Required  Quantity 

Units 

Latitude  of  point  "A"  (takeoff) 

decimal 

deg 

Longitude  of  point  "A"  (takeoff) 

decimal 

deo 

Latitude  of  point  "B"  (landing) 

decimal 

deg 

Longitude  of  point  "B"  (landing) 

decimal 

deg 

Julian  base  date  (1-366) 

none 

Greenwich  base  time  (0.0-24.0) 

hr 

Forecast  hours  ahead 

hr 

Altitude  index  (1  =  25,000  ft,  2  =  35,000  ft) 

none 

Wind  option  or  mode  (0  *  calm  wind. 

2  =  90-percent  worst  wind,  3  *  mean  wind) 

none 

Airspeed 

kt 

2.4  Restrictions  Imposed 

The  WFSM  was  originally  constrained  to  operate  within  1.5K  words  of  core  storaqe. 
This  restriction  was  not  met.  The  model  now  requires  5. IK  plus  file  control  block 
and  buffer. 

2.5  Future  Requirements 

Simulation  models  such  as  MAC'S  COLOSSUS, in  which  the  WFSM  must  reside,  have  to 
approximate  the  major  actions  and  decisions  made  in  operational  practice.  Some  of 
these  operator  decisions  are  based  on  climatology,  some  on  forecasts,  and  some  on 
observations.  In  the  present  model  the  day-to-day  wind  "observations"  are  always 
mean  or  cl imatologically  expected  values.  There  is  no  day-to-day  variability  in  the 
"observed"  winds.  Furthermore,  the  model  cannot  simulate  wind  or  wind  factor 
"forecasts. " 

When  the  "observed"  winds  always  correspond  to  a  climatologically  expected 
value  the  upper  atmosphere  wind  flow  does  not  change.  It  is  more  realistic  to  permit 
a  dav-to-day  variability  in  the  wind  factor,  a  variability  that  is  the  statistical 
analog  of  the  natural  evolution  of  upper-air  wind  patterns.  Future  requests  will 
doubtless  state  a  requirement  for  "simulated"  wind  factors  having  a  natural 
variability.  The  architecture  of  subroutine  ENRWND  will  accomodate  a  "simulated" 
wind  factor  technique.  The  present  subroutine  defaults  to  the  mean  wind  factor  case. 


Future  requests  may  include  a  requirement  for  "forecast"  winds.  If  climatology, 
"forecasts,"  and  "observations"  are  all  realistically  included  in  the  WFSM,  users  of 
that  model  can  determine,  in  the  larger  context  of  their  own  mission  simulations,  the 
optimum  mix  of  climatology,  forecasts,  and  observations. 


Chapter  3 


WIND  FACTOR  SIMULATION  MODEL 


3.1  General  Description  of  the  Method  and  Overall  Simulation  Logic 


Sawyer's  eauivalent  headwind  technique  (Air  Weather  Service,  1977)  is  combined 
with  the  mathematics  of  qreat  circle  navigation  in  a  constellation  of  seven  FORTRAN 
subprograms  whose  purpose  is  given  in  Table  3. 


TABLE  3.  WIND  FACTOR  SIMULATION  MODEL  SUBPROGRAMS. 


Subprogram 

ENRWND 


DISTAN 

SPHGLO 

GRTCIP 

HDG 

BRLNG 

BRLAT 


Purpose 


Main  enroute  wind  factor  subprogram  called  directly  by  user  to 
compute  ground  speed  from  given  airspeed  in  anv  of  three  modes 
for  two  flight  levels  and  four  seasons  of  the  year. 

Called  by  ENRWND  and  HDG  or  directly  bv  user  to  comoute  qreat 
circle  distance  in  NM  between  any  two  points  on  the  qlobe. 

Conversion  of  spherical  coordinates  to  global  latitude  and  longi¬ 
tude  or  vice  versa.  Called  by  GRTCIR. 

Solution  of  great  circle  eauation  for  latitude  given  longitude, 
or  for  longitude  given  latitude.  Called  by  ENRWND. 

Calculation  of  initial  heading  along  a  qreat  circle  route  flown 
from  a  given  oriqin  to  a  given  destination.  Called  by  ENRWND 
and  GRTCIR. 

For  a  spherical  grid  system  whose  longitude  boundaries  are  spaced 
at  30°  intervals,  finds  the  longitude  grid  values  bracketing  a 
given  arbitrary  lonaitude.  Called  by  ENRWND. 

For  a  spherical  grid  system  whose  latitude  boundaries  are  spaced 
at  15°  intervals,  finds  the  latitude  qrid  values  bracketing  a 
given  arbitrary  latitude.  Called  by  ENRWND. 


The  keystone  subprogram  of  the  Wind  Factor  Simulation  Model  is  ENRWND,  which 
handles  most  of  the  aircraft  navigation  except  solution  of  the  eouation  of  a  qreat 
circle  and  also  accomplishes  all  of  the  meteorological  computations  of  the  wind 
factor. 

The  ENRWND  subprogram  naviqates  a  simulated  aircraft  over  a  great  circle  route 
from  point  "A"  to  point  "B"  in  a  global  (latitude/longitude  oriented)  a  rid  system 
whose  resolution  is  15°  latitude  by  30°  longitude.  The  subroutine  determines  the 
grid  sector  location  of  the  aircraft  and  then  determines  the  length  of  the  aircraft’s 
oath  through  the  sector.  Next,  the  routine  consults  the  wind  data  base  for  the 
particular  grid  sector,  flight  level,  and  season  to  obtain  sector-mean  wind 
direction,  speed,  and  variance.  Cross-track  and  along-track  components  of  the  wind 
are  computed.  These  components  are  weighted  by  the  length  of  the  aircraft's  track 
through  the  sector  and  are  accumulated.  At  the  end  of  the  simulated  flight,  the 
accumulators  are  divided  by  the  total  GCD  for  the  flight,  producing  a 
distance-weighted,  route-mean  cross-track  wind  component,  an  along-track  component, 
and  a  variance.  From  these  components  a  route-mean  wind  factor  is  computed.  If  the 
90-percent  worst  wind  option  has  been  selected,  that  wind  factor  is  statistically 
adjusted  to  the  10-percent  risk  value.  The  wind  factor  is  then  added  to  the 
airspeed  to  obtain  the  route-mean  airspeed.  If  the  calm  wind  option  is  requested, 
all  of  these  computations  are  bypassed,  and  a  zero  wind  factor  is  used.  Reauests  for 
the  simulated  wind  option  default  to  the  mean  wind.  In  this  case,  the  mean  wind  is 
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developed  directly  from  the  mean  wind  data  base.  In  the  90-percent  worst  case,  this 
mean  wind  factor  is  statistically  adjusted  to  a  value  such  that  only  10  percent  of 
fliqhts  over  that  route  will  experience  a  worse  wind  factor  (10-percent  risk). 

3.2  Sawyer's  Equivalent  Headwind  Technique 

Leq  and  overall  route  wind  factors  ate  computed  by  Sawver's  equivalent  headwind 
technique  (as  described  in  Section  3.1). 

A  wind  factor  W  can  be  defined,  as  the  difference  between  the  qround  speed  |g!  of 
an  aircraft  and  its  true  airspeed  |a  |  •  This  relationship  is  qiven  in  the  equation: 

W  =  \t\  -  IaI  =  G  -  A  (1) 

Sawyer  found  that  a  mean  or  climatoloqical  overall  wind  factor  for  a  route  is 
qiven  by  the  relationship: 

[w]  =  [u]  -  sj  (tv2]  +  ^r1)  {2) 

where  []  indicates  an  average  over  all  the  legs,  and  where 
W  =  mean  wind  factor 

U  =  mean  wind  component  along  aircraft  track 
V  =  mean  crosswind  component 

• 

o2  =  mean  variance 

The  first  term  to  the  riqht  of  the  equal  sign  in  Eo.  (2)  is  the  tailwind  component, 
and  the  second  term  is  the  effect  of  the  crosswind. 

In  calculating  the  90-percent  worst  wind  factor  (i.e.,  10-percent  risk),  Sawver 
used  the  equation: 

W90#  =  ^  "  1,28  °w  (3> 

where  ow  is  related  to  C cr J  by  the  Sawyer  K-factor.  This  K-factor  is  route  lenqth 
dependent  and  is  expressed  as 


K 


(4) 


Sawyer  provides  a  table  of  K  as  a  function  of  route  length.  For  convenience,  a 
geometric  curve  was  fitted  to  Sawyer's  tabular  data  using  the  Hewlett-Packard  65  Stat 
Pac  2  curve  fit  program  (Hewlett-Packard  Company,  1975c).  The  curve  is 


K  **  0.70585066  x  (0. 99983056 )d 


where  d  *  great  circle  distance  in  nautical  miles. 

3.3  Grid  System 

The  VFSH  must  be  economical  in  terms  of  computer  core  storage  and  run  time: 
yet  it  must  provide  a  global  wind  factor  capability  to  the  user.  These  conflictinq 
reguirements  dictated  the  use  of  a  coarse  grid.  The  original  idea,  motivated  bv  the 
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^*1 


latitude/lonqitude  orientation  of  the  winds  used  as  input  data,  was  to  have  a  grid 
system  in  global  (latitude/lonqitude)  coordinates.  A  qrid  would  have  a  spacing  of 
15°  latitude  by  30°  longitude  and  would  cover  the  qlobe,  reauirinq  144  qrid  sectors 
(12  sectors  from  pole  to  pole  times  12  sectors  around  the  eauator).  Closet 
investigation  of  the  latitudes  of  the  COLOSSUS  terminals  and  routes  indicated  that 
qrid  sectors  were  not  necessary  above  75°N  and  below  60°S.  Since  three  quantities 
(wind  direction,  speed,  and  variance)  are  stored  for  each  grid  sector  for  four 
seasons  and  two  flight  levels,  approximately  0.87K  words  of  core  storaqe  could  be 
saved  by  reducing  the  number  of  qrid  sectors  from  144  to  108. 

The  grid  system  chosen  for  implementation  is  in  global  (latitude/lonqitude) 
coordinates,  with  a  resolution  of  15°  latitude  by  30°  longitude.  The  qrid  system 
extends  from  75°N  to  60°S  (9  rows)  and  qlobally  in  longitude  (12  columns).  Hence, 
there  are  108  grid  sectors  as  shown  in  Figure  1.  Note  that  column  1  is  keyed  on  30  W 
not  0°. 
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Figure  1.  Grid  System. 

Because  this  grid  system  does  not  include  the  Doles,  the  Wind  Factor  Simulation 
model  cannot  consider  over-the-pole  flights.  The  coarseness  of  the  qrid  system 
orecludes  use  of  this  model  for  calculation  of  operational  wind  factors.  The  coarse 
orid  is  suitable  for  simulation,  however,  and  could  be  made  finer  if  improved 
meteorolooical  resolution  is  desired. 

3.4  Data  Base 


Basic  data  for  the  WFSM  are  the  USAFETAC  Sinqle  Xnteorated  Operations  Plan  (STOP) 
winds.  With  a  period  of  record  extending  from  January  1972  to  December  1976,  the 
SIOP  winds  contain  mean  u-component,  v-component,  and  vector  standard  deviation 
tabulated  on  a  5-deqree  offset  qrid.  This  qrid  is  much  finer  than  that  used  by  the 
WFSM. 

To  prepare  a  data  base  for  the  Wind  Factor  Simulation  Model's  15°  latitude  by  30° 
longitude  grid  system,  simple  averaging  was  used.  The  u-component,  the  v-component, 
and  the  vector  standard  deviation  were  averaged  separately.  The  average  was 
performed  by  extracting  from  the  SIOP  winds  three  values  latitudinally  and  six  values 
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longitudinally.  The  resulting  3  x  6  *  18  values  were  summed  and  the  sum  divided  by 
18  to  obtain  the  mean  data  values  for  each  grid  sector.  In  this  way,  mean 
u-components,  v-components,  and  vector  standard  deviations  were  obtained  for  all  108 
grid  sectors.  Such  a  data  base  was  constructed  for  January  (winter),  April  (spring), 
July  (summer),  and  October  (fall)  for  altitudes  25,000  feet  (taken  from  400-mb  winds) 
and  35,000  feet  (taken  from  250-mb  winds).  Using  an  in-house  USAFFTAC  program  WIND, 
the 


4  season  x  2  altitude  =  8  sets 

of  108-sector  wind  data  were  converted  from  u-components,  v-components,  and  vector 
standard  deviation  to  direction  (beta  angle,  described  in  Appendix  B)  ,  speed,  and 
variance.  After  conversion,  the  winds  were  stored  seauentially  in  a  data  file  in  the 
order  described  in  Table  4.  The  data  file  consists  of  864  records  (4  seasons  x  2 
altitudes  x  108  sectors).  Each  sector  of  the  data  file  requires  one  record  of  disk 
storage  containing  the  three  elements  of  information  shown  in  Table  5. 


TABLE  4.  STRUCTURE  OF  WIND  DATA  FILE. 

Flight  Level 

Number  of 

Record 

Season 

(hundreds  of  feet) 

Sectors 

Numbers 

Winter 

250 

108 

1 

-  108 

350 

108 

109 

-  216 

Spr inq 

250 

108 

217 

-  324 

350 

108 

325 

-  432 

Summer 

250 

108 

433 

-  540 

350 

108 

541 

-  648 

Fall 

250 

108 

649 

-  756 

350 

108 

757 

-  864 

TABLE 

5.  FORMAT 

OF  WIND  DATA  RECORD 

FORTRAN 

ELEMENT 

UNITS 

COLUMNS 

FORMAT 

Wind  Direction 
(Beta  Angle) 

radians 

1-15 

F15.5 

Wind  Speed 

knots 

16  -  30 

F15.5 

Variance 

knots2 

31  -  45 

F15.5 

To  be  used  by  the  Wind  Factor  Simulation  Model,  the  data  base  must  be  loaded  from 
the  data  file  into  the  arrays: 

DIP (870)  Wind  direction  beta  angle  (radians) 

SPD(870)  Wind  speed  (knots) 

VAR (870)  Wind  variance  (knots2) 

These  arrays  are  located  in  the  labeled  COMMON  block  WEA  and  require  2.6K  words  of 
computer  core  storage.  The  COMMON  block  must  be  loaded  with  wind  data  before  the 
WFSM  is  first  executed.  Usually  an  initializing  routine  external  to  the  WFSM  is  used 
for  this  purpose.  That  routine  reads  864  wind  records  in  order  and  stores  them, 
without  rearrangement,  in  the  COMMON  block  WEA.  Thereafter,  the  wind  data  base  can 
be  accessed  in  terms  of  a  record  number  NRND  as  follows: 

DIR(NRND)  SPD(NRND)  VAR (NRND) 
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The  record  number  NRND  is  computed  by: 

NRND  =  (Season  Index  -  1)  *  (Total  Number  Altitudes) 

*  (Total  Number  Grid  Sectors) 

+  (Altitude  Index  -  1)  *  (Total  Number  Grid  Sectors) 

+  (Sector  Number)  (5) 


3.5  Aircraft  Navigation  Technique  and  Logic 

Aircraft  navigation  is  handled  principally  by  the  subroutine  ENPWND.  The 
subroutine  subprogram  GRTCIR  solves  the  equation  of  a  great  circle. 

The  simulated  aircraft  is  navigated  along  a  great  circle  route  between  points  "A" 
and  "B"  by  moving  from  a  point  called  the  "Current-Point"  (XLATC,  XLNGC)  to  a  point 
called  the  "Next-Point"  (XLATNP,  XLNGNP) .  Most  of  the  computation  of  the  Next-Point 
is  done  in  global  (latitude/longitude)  coordinates. 

For  any  given  Current-Point,  the  Next-Point  is  either  the  point  where  the  great 
circle  route  crosses  a  latitude  grid  line  (0,  15,  30,  ...,  deg)  or  the  point  where  a 
route  crosses  a  longitude  grid  line  (0,  30,  60,  ...,  deg),  whichever  is  closer  to  the 
Current-Point.  This  is  illustrated  in  Figure  2. 


Figure  2.  Aircraft  Navigation  Technique. 
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At  the  Current-Point  the  aircraft's  heading  can  he  calculated  via  the  headina 
formula  (function  HDG) .  Knowledge  of  the  location  of  the  Current-Point  and  the 
aircraft's  heading  at  that  point  dictates  the  next  latitude  crossino  and  the  next 
longitude  crossing.  Candidate  #1  for  the  Next-Point  is  the  point  of  next  latitude 
crossing,  XLATNX  (and  associated  longitude,  XLNGAS) .  Candidate  #2  is  the  point  of 
next  longitude  crossing,  XLNGNX  (and  associated  latitude,  XLATAS) . 

Calculation  of  lonqitude,  XLNGAS,  associated  with  a  oiven  latitude,  XLATNX,  is 
done  iteratively  by  subroutine  GRTCIR,  which  solves  the  equation  of  a  qreat  circle. 
Newton's  iterative  method  with  an  exact  spherical  triangle  first  quess  is  used. 
Virtually  all  the  complexity  of  subroutines  GRTCIR  and  ENRWND  is  attributable  to 
solvinq  the  eauation  of  a  great  circle  for  longitude,  qiven  the  latitude. 

Calculation  of  latitude,  XLATAS,  associated  with  a  given  lonqitude,  XLNGNX,  is 
accomplished  by  subroutine  GRTCIR. 

Candidates  for  Next-Point  are  then  as  follows 

Candidate  #1  (XLATNX,  XLNGAS) 

Candidate  «2  (XLATAS,  XLNGNX) 

With  Candidates  #1  and  |2  identified,  subroutine  DISTAN  is  then  used  to  calculate 
great  circle  distances  as  follows: 


SI 

GCD 

from 

Current-Point 

to 

Candidate 

#1 

S2 

GCD 

from 

Current-Point 

to 

Candidate 

#2 

The  candidate  with  the  shorter  distance,  SI  or  S2,  is  selected  as  the  Next-Point. 

The  Current-Point  can  then  be  thought  of  as  the  point  where  the  aircraft  enters 
the  grid  sector,  and  the  Next-Point  as  the  point  where  the  aircaft  departs  the 
sector.  The  distance  between  Current-Point  and  Next-Point,  SNP,  is  the  distance  by 
which  the  wind  components  and  variance  are  weighted  in  computing  the  route-mean  wind 
factor . 


3.6  Meteorological  Technique  and  Logic 

Subroutine  ENRWND  calulates  the  meteorological  ouantities  reouired  by  the  Wind 
Factor  Simulation  Model. 

With  Current-Point  and  Next-Point  established,  an  appropriate  midpoint  (XLATM, 
XLNGM)  between  them  is  computed.  The  number,  NSEC,  indicating  which  of  the  108  grid 
sectors  through  which  the  aircraft  is  passing  is  then  calculated  based  on  (XLATM, 
XLNGM).  See  the  description  of  the  grid  system  in  section  3.3  above.  Since  four 
seasons  and  two  flight  levels,  as  well  as  108  grid  sectors  are  provided  for,  a  datum 
number  (NRND)  is  obtained  from  the  grid  sector  number  (NSEC,  values  1-108),  the 
altitude  index  (IALT,  values  1-2)  and  the  season  index  (ISEASN,  values  1-4)  using  Ea. 
(6)  . 

With  the  datum  number,  NRND,  calculated,  the  data  base  (see  description  in 
section  3.4)  is  consulted  to  obtain  sector-mean  wind  direction  (DIRM,  radians  in  beta 
angle  form,  as  described  in  Appendix  8),  speed  (SPDM,  knots),  and  variance  (VARM, 
knots2) : 


DIRM  -  DIR (NRND) 
SPDM  ■  SPD(NRND) 
VARM  -  VAR (NRND) 
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Function  HDG  computes  the  aicctaft  heading  (HDGM)  at  (XLATM,  XLNGM)  in  beta  anale 
format  {see  Appendix  B) .  In  that  format  HDGM  is  represented  by  the  angle  ALPHA. 

ALPHA,  the  wind  direction,  DIRM,  and  the  wind  speed,  SPDM,  are  used  to  calculate 
the  along-ground-track  wind  component,  VGM,  and  the  cross-ground-track  wind  compo¬ 
nent,  VCM: 


VGM  -  SPDM  X  COS  (ABS (DIRM  -  ALPHA)) 


(5) 


VCM  -  SPDM  X  SIN  (DIRM  -  ALPHA) 


(9) 


A  time-averaged  wind  factor  for  the  present  leg  (WPAR)  is  then  computed  usina 
Sawyer's  equivalent  headwind  formula.  Eg.  (2).  No  further  use  is  made  of  WBAR. 

The  quantities  VGM,  VCM,  and  VARM  are  multiplied  by  the  sector  or  leg  lenqth, 
SNP,  and  the  product  is  accumulated  in  VGACC,  VCACC,  and  VACC.  Then  the  Next-Point 
becomes  the  Current-Point,  and  a  new  Next-Point  is  calculated.  The  process  is 
repeated  for  the  next  leg  of  the  simulated  flight. 

After  the  last  leg  of  the  flight  has  been  processed,  the  accumulators  VGACC, 
VCACC,  and  VACC  are  divided  by  the  total  route  length.  The  resultant  distance- 
weighted  along-track,  cross-track,  and  variance  quantities  are  then  used  in  the 
Sawyer  equivalent  headwind  formula,  Eq.  (2),  to  obtain  the  distance-weighted, 
route-mean  wind  factor,  WBARBR.  To  obtain  the  90-percent  worst  wind  factor,  WBARBR 
is  statistically  adjusted  to  the  90-percent  worst  value  by  Sawyer's  technique. 

As  part  of  the  wrap-up  actions  of  subroutine  ENRWND,  the  route-mean  wind  factor, 
WBARBR,  is  added  to  the  airspeed,  ASPEED,  to  get  a  ground  speed,  GSPEED,  for  return 
to  the  program  calling  subroutine  ENRWND. 


3.7  Equation  of  a  Great  Circle 


In  the  wind  factor  simulator,  great  circle  routes  are  flown  from  point  "A,"  the 
origin,  to  point  "B,"  the  destination.  A  great  circle  is  the  intersection  of  a 
sphere  and  a  plane  passing  through  the  center  of  the  sphere,  shown  in  Figure  3.  As  a 
mathematical  curve,  the  great  circle  must  be  representable  in  terms  of  an  equation. 
One  such  equation  may  be  developed  from  vector  algebra  as  follows. 

Consider  that  is  the  position  vector  from  the  center  of  the  sphere  to  point 
"A"  along  the  great  circle,  while  QB  is  the  position  vector  of  point  "B."  The  plane 
of  the  great  circle  is  defined  by  the  vector  cross  product: 


where  (Jr  is  normal  to  the  plane. 

Let  6  be  the  position  vector  of  any  intermediate  point  "D"  along  the  great 
circle.  Since  0?  is  normal  to  0,  an  equation  describing  the  great  circle  is 


5  •  ($A  *  V  "  0 


(11) 


Points  "A",  "B",  and  "D"  are  ordinarily  described  in  terms  of  latitude  md 
longitude,  which  are  here  celled  global  coordinates.  Global  coordinates  are  easily 
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and  where 


The  cross  product  is 


QBx  -  r  sin  0B  cos  ^  -  rqfix 

Qgy  -  r  sin  eB  Sin  ^  -  rq^ 

<3bz  -  r  cos  0B  -  rqBz 

Qx  m  r  sin  0  cos  «p  ■  rqx 

Qy  *  r  sin  0  sin  «p  ■  rq^ 

-  r  cos  0  -  rq 


qBx  *  sin  ®b  008  (*IB 


qBy  "  8in  eB  8in  *B 


cos  0_ 


qx  ■  sin  0  cos  <p 


qy  ■  sin  0  cos 


q2  ■  cos  0 


-  $A  x  5b  -  rqAx 


t 

1 

k 

^Ax 

rqAy 

rqAz 

rqBx 

rqBy 

rqBz 

qAx  “ 

sin 

®A 

cos 

*A 

qAy  * 

sin 

0A 

sin 

*A 

1 

,  (16) 

qAz  " 

cos 

«A 

J 
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or 


[^Ay  qBz  "  qAz  qBy}  *  +  (qAz  qBx  “  qAx  qBz>  ^ 
+  (qAx  qBy  "  qAy  qBx^ 

5ip  “  r  (^Tx  ^  +  qTy  ^  +  qTz  k^ 

where 

qTx  ”  qAy  qBz  "  qAz  qBy 
qlV  “  qAz  qBx  "  qAx  qBz 
qTz  *  qAx  qBy  ”  qAy  qBx 
The  vector  eauation  of  the  great  circle  is 

Q  •  «  0 

or  in  component  form, 

qx  <hx  +  qTy  +  qz  qTz  "  0 


(20) 


(21) 


(22) 


(23) 


where  g>rx,  <*Ty»  and  ®Tz  are  known  from  points  "A"  and  "B."  Substituting  for  gx,  Qy 
and  az  qives  the  results, 

/  n  j  1 

(sin  8  cos  cp)  +  (sin  8  sin  <p)  q^y  +  (cos  8)  QTz  ”  0 

which  is  the  equation  of  a  great  circle  in  rectangular  component  form. 

Two  solutions  of  the  eauation  of  a  great  circle  are  possible: 

Case  1:  Solve  for  9  with  cp  known 
Case  2:  Solve  for  cp  with  6  known 

For  Case  1,  a  deterministic  solution  can  be  obtained  by  dividing  the  great  circle 
equation  by  sin  8  • 


i 

J 


/ 


or 


cos  ®  +  q^y  sin  cp  +  ril'  -s  *  0 


tan  l" 


0  “  tan  ^tx  008  *  +  ^  8in  ^ 


(25) 


(26) 
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which  is  easily  evaluated  and  solved,  since  q>px,  aT v»  QTz»  and  <P  are  known. 

Case  2  proves  more  difficult  to  solve  in  that  it  takes  the  form  of  a  transcen¬ 
dental  eauation.  Newton's  iterative  method  (sometimes  called  the  Newton-Raphson 
method)  for  the  solution  of  nonlinear  algebraic  equations  was  selected  as  the 
numerical  technique  with  which  to  solve  the  great  circle  eauation  for  tp  ,  where 
6  is  known. 

3.8  Newton's  Iterative  Method 

The  zeros  or  roots  of  the  nonlinear  function, 

f(g>)  =  0  <27‘ 

in  Newton's  method,  are  found  by  the  recursion  relation, 

„("«)  .  „<»> .  jasjju  ,2., 


where  n  is  the  iteration  count  and  f'  is  the  first  derivative  of  the  function.  In 
the  case  of  the  great  circle  formula  in  the  form  of  Eg.  (27), 

f(<p)  -  q-Tx  cos  9  +  Qry  sin  <P  +  (29: 


f '  M  ■  "qTx  sln  9  +  qTy  006  9 


Newton's  method  requires  a  first  guess,  cp(0).  The  first  guess  must  be  "close 
enough"  to  the  real  root  to  insure  convergence  of  the  method,  which  then  is  quite 
rapid,  with  the  number  of  good  decimal  digits  nearly  doubling  at  each  iteration.  If 
the  first  guess  is  not  close  enough,  Newton's  method  may  diverge.  For  the  great 
circle  equation,  the  first  guess,  g>(0) ,  must  be  within  about  35°  (0.61  radians)  of 
the  true  root  to  insure  convergence.  Methods  for  generating  the  first  guess  are 
discussed  in  the  following  section  of  this  report. 

3 . 9  First  Guess  Methods  for  Newton's  Iterative  Technique 


Solving  the  transcendental  equation  for  longitude  (azimuth),  given  the  latitude 
(colatitude),  by  Newton's  iterative  technique  reauires  a  first  guess  for  the 
azimuth  <P  .  In  practice,  obtaining  this  first  guess  is  critical  for  convergence 
of  “ewton's  technique,  which  otherwise  may  diverge  or  converge  to  false  roots. 

In  subroutine  GRTCIR,  three  methods  are  potentially  used  to  obtain  the  first 
guess  azimuth  tp9.  The  first  such  method,  not  normally  used  except  in  debuoaing,  is 
for  the  user  to  specify  the  value  of  tp9  .  The  second  method  is  that  of  approximately 
planar  trianqles.  This  method  has  important  limits  on  its  reliability  and  so  is  used 
only  to  find  the  longitude(s)  where  the  great  circle  crosses  the  equator.  The  third, 
and  most  exact  first  guess  method,  is  that  of  the  exact  spherical  triangle.  The 
spherical  triangle  method  requires  knowledge  of  the  longitude  of  equator  crossing  and 
thus  cannot  be  used  initially. 

In  application  for  any  given  route,  the  approximate  planar  triangle  method  is 
used  initially  to  obtain  the  longitude  of  equator  crossing.  Thereafter,  the  exact 
spherical  triangle  method  is  used  throughout  the  route. 

3.9.1  Method  of  Approximately  Planar  Triangle.  A  segment  of  a  great  circle  is  shown 
by  the  line  PF  In  Figure  4,  where  the  point  E  may  be  either  the  origin  A  of  the 
flight  or  the  destination  B,  whichever  is  closer  to  the  equator.  If  E  coincides  with 
A,  then  F  must  coincide  with  B?  if  E  coincides  with  B,  then  F  must  coincide  with  A. 
If  the  route  is  extended  to  the  equator,  the  equator  crossing  point  A'  is  defined  as 
shown  in  Figure  4. 
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Figure  4.  Approximately  Planar  Triangle. 


The  spherical  triangle  A'EE'  must  be  solved  for  the  arc  A 'E  ■  in  order  to  obtain 
the  azimuth  ®A>  of  the  equator  crossing  point  A'.  Not  enough  infomation  is  known 
about  this  trianqle  to  solve  it  spherically.  If  A  'EE '  is  approximated  by  a  planar 
right  triangle,  an  approximation  to  the  distance  A'E’  can  be  obtained  through  the 
tanoent  of  the  angle  E,  which  is  known  from  the  heading  of  the  aircraft  at  point  E: 


A'E'  «*  E>'fe  tan  E  (radians)  (31) 


Put  the  angular  distance  E  *E  is  simply  the  known  latitude  of  point  E  in  radians. 
Hence: 


nrr  "  LATj.  tan  E  (radians)  (32) 


and  the  first  guess  eauator  crossing  azimuth  is 


-  v  -  r*r  (radians)  (33) 


The  first  guess  is  then  improved  to  the  final  value  of  by  Newton's 
iterative  technique,  described  in  paragraph  3.8  above. 
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3*5*2  Method  of  Exact  Spherical  Triangle.  A  segment  of  a  great  circle  is  shown  as 
the  line  AB  of  Figure  5,  where  point  A  is  the  origin  of  the  flight,  and  point  B  is 
its  destination.  Point  D  is  any  intermediate  point  along  the  AB  great  circle.  The 
riqht  spherical  triangle  of  interest  is  A'DE' .  where  point  A'  is  the  eauator  crossing 
established  above,  and  the  angular  distance  DE 1  is  given  by  the  known  latitude  of  the 
intermediate  point  D. 


Figure  5.  Spherical  Triangle  Opening  Eastward. 


Since  the  position  of  A'  is  known,  the  angle  A’  can  be  obtained  from  the  heading 
of  the  aircraft  at  A'.  The  unknown  of  interest  is  the  great  circle  angular  distance 
A’D. 

The  law  of  sines  for  spherical  triangles  gives  the  result: 


(34) 


(radians)  (35) 


sin  E'  m  sin  A 
sin  A' D  sin  DE7 


rrr  .  -1  /sin  E'  •  sin  DE' 

A  D  ’  sln  V - sinT - > 


From  the  spherical  right  triangle  it  is  known  that  the  angle  E*  is  90°:  hence: 


-rrr-  .  -1  /sin  DE'  > 
A  D " 8ln 


(radians) 
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for  which  there  are  solutions  in  all  quadrants.  In  practice,  the  principal  solution 


TT 

•  ? 


<  A'D, 


(37) 


and  the  solution, 


(radians) 


(38) 


are  of  interest. 

It  is  desirable  to  estimate  the  angular  distance  £TET.  Since  the  azimuth  <Fa'  of 
the  equator  crossing  point  A'  is  known,  the  problem  of  obtaining  A 1 ET  reduces  to  that 
of  estimating  the  azimuth  epp  of  the  point  D,  which  is  the  same  as  the  azimuth  <Pe' 
of  the  point  E*.  The  problem  can  be  addressed  in  terms  of  longitude  via  the  great 
circle  distance  formula: 


cos  (A* D)  -  sin  (LATA» )  sin  (LATD) 

cos  (LNGd  -  lngaO  -  - cos  (UTATT ~os'TD^ -  <39) 

But  the  latitude  of  the  equator  is  0°,  so 

LNGd  -  LNOA'  ±  cos"1  [cos  (49) 

where  the  negative  sign  is  used  for  the  spherical  triangle  opening  eastward  (Figure 
5),  and  the  positive  sign  is  used  for  the  spherical  triangle  ooening  westward  (Fiaure 


Figure  6.  Spherical  Triangle  Opening  Westward. 
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Figure  7.  Route  Crossing  Equator  Shows  Reversal  of  Spherical  Triangle. 


The  criterion  for  "ooening  eastward"  is  that  point  S'  must  be  east  of  ooint  A*. 
In  a  triangle  that  "opens  westward",  ooint  E'  is  west  of  ooint  A'.  Routes  that  cross 
the  eouator  exhibit  reversal  of  the  spherical  triangle  direction  at  the  point  of 
equator  crossing,  as  shown  in  Figure  7. 

3.10  Great  Circle  Distance 

The  Wind  Factor  Simulation  Model  (WFSM)  computes  the  great  circle  distance  (GCD) 
through  its  subroutine  DISTAN.  DISTAN  was  developed  from  a  program  of  the  Hewlett- 
Packard  Company  (HP).  The  HP-65  Navigation  Pac  1  (Hewlett-Packard  Company,  1975b) 
describes  the  method. 

The  basic  equation  for  the  great  circle  distance  in  radians  is 


GCDr  ■  cos-1tsin(LATA)sin(IATB)+cos(LATA)cos(LATB)cos(LNGB-LNGA)]  (41 


where 

LATA,  LNGa  «  coordinates  of  initial  point  "A"  in  radians 
tiATg,  LNGjj  ■  coordinates  of  final  point  "B"  in  radians 

GCDr  ■  great  circle  distance  in  radians  from  "A"  to  "B" 


To  convert  to  GCD  in  nautical  miles  (NM) ,  the  relation 


GCD_  =  a  •  GCD 
o  r 


(42) 


is  used,  where  "a"  is  the  radius  of  the  earth  in  nautical  miles  (approximately  3440 
NM)  . 

DISTAN  is  employed  for  more  than  just  computing  GCDs  from  departure  to 
destination.  It  is  called  by  subroutine  ENRWND  every  time  the  aircraft  is  ready  to 
move  from  its  Current-Point  to  the  Next-Point  along  the  great  circle  track.  The 
Next-Point  is  always  the  candidate  point  with  the  shortest  GCD  between  it  and  the 
Current-Point.  DISTAN  can  be  called  as  few  as  four  times  (route  does  not  cross  qrid 
boundaries)  or  as  many  as  25  or  more  times  depending  upon  route  length  and  number  of 
grid  crossings. 

3.11  Great  Circle  Heading 


The  WFSM  computes  the  initial  heading  (degrees  clockwise  from  north  toward  which 
the  aircraft  is  moving)  along  a  great  circle  route  from  point  "A"  to  point  "B"  by 
means  of  the  function  subprogram  HDG.  This  function  was  developed  from  the 
Hewlett-Packard  65  Aviation  Pac  1,  Great  Circle  Navigation  (Hewlett-Packard  Company, 
1975a)  .  The  function  uses  input  latitudes  and  longitudes  for  points  "A"  and  "B"  and 
converts  them  from  degrees  to  radians,  then  uses  the  heading  formula. 


_ -i  rsin  (“V  -  sin  (“V  cos  (GCDrh 

HDGr  =  cos  [ - sin'  (Gfi£>r)  cos  (IMaJ - J 


(43) 


to  find  the  initial  heading,  HDGr ,  in  radians  at  point  "A".  The  other  variables  are 
described  in  paragraph  3.10  above.  The  great  circle  distance  in  radians,  GCDr ,  is 
computed  by  use  of  subroutine  DISTAN,  which  gives  GCD„,  in  nautical  miles.  GCDm  is 
converted  to  GCDr  by  the  relation. 


GCDr  = 


gcd„, 

m 

a 


(44) 


where  "a"  is  the  radius  of  the  earth  in  nautical  miles  (approximately  3440  NM) . 
If 


sin  (LNGa  -  LNGb)  <0  (45) 


then  the  first  auadrant  arc  cosine  solution  is  not  desired.  A  correction  is  made  by 
sett  ing 


HDG'  =  2tt  -  HDG  (46) 

r  r 

Then  the  heading  in  radians  is  converted  to  the  heading  in  degrees  by  multiplying  by 

180/  TT  . 
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Chapter  4 


MODEL  CAPABILITIES  AND  LIMITATIONS 


4 . 1  Model  Capabilities 


The  Wind  Factor  Simulation  Model  (WFSM)  is  capable  of  producina  mean  overall  wind 
factors  for  great  circle  routes.  It  does  so  in  any  of  three  modes:  calm  wind  case, 
90-percent  worst  wind  case,  and  the  mean  wind  case.  The  wind  factor  can  be  produced 
for  two  flight  levels  and  for  the  four  seasons.  In  addition  the  model  can,  throuqh 
its  subroutine  DISTAN,  provide  the  qreat  circle  distance  (GCD)  between  any  two  points 
on  the  qlobe.  Since  it  produces  both  wind  factors  and  GCD,  it  is  capable  of  calcu¬ 
lating  both  the  ground  speed  and  the  adjusted  flying  time  between  points  "A"  and  "D" 
if  given  the  airspeed. 

4 . 2  Model  Assumptions 


The  WFSM  assumes  that  all  routes  are  great  circle  routes  or  can  be  decomposed 
into  several  legs  each  of  which  is  a  great  circle.  Further,  it  assumes  that  climb 
winds  and  descent  winds  play  a  negligible  role  in  determinina  the  route-mean  wind 
factor;  the  aircraft  is  always  "at  altitude".  Since  an  arithmetic  rather  than  a 
harmonic  mean  is  used  in  computing  the  route-mean  wind  factor,  it  is  assumed  that  the 
wind  speed  is  less  than  or  equal  to  one-third  of  the  airspeed. 

4.3  Model  Limitations 


The  WFSM  cannot  produce  a  "forecast"  wind  factor  or  a  "simulated"  wind  havinq  a 
day-to-day  variability.  Bence,  within  any  qiven  season  and  for  any  given  altitude,  a 
particular  route  will  alwavs  experience  the  same  wind  factor,  reqardless  of  the 
passage  of  time.  Mean  winds  are  always  substituted  for  simulated  and  "forecast" 
winds.  Temperature,  aircraft  performance,  fuel  consumption,  and  other  factors 
considered  by  typical  flight  planning  models  are  not  included  here  since  the  present 
model  deals  with  wind  only. 

Operations  north  of  75°N  or  south  of  60°S  are  not  permitted. 

The  coarseness  of  the  present  15°-latitude  by  30°-longitude  grid  system  precludes 
use  of  the  model  for  calculation  of  operationally  realistic  wind  factors.  Some 
simple  modifications  to  the  existing  data  base  and  grid  system,  along  with 
corresponding  chanqes  in  the  software,  can  remove  this  limitation. 

The  followina  restrictions  must  be  adhered  to: 


a.  No  route  should  be  flown  directly  over  either  pole. 

b.  Neither  pole  should  serve  as  a  point  "A"  or  "B." 


c . 

be  aborted. 


Circumferential  or  round  robin  flights  in  which  "A"  and  "B"  coincide  will 
Break  such  flights  into  smaller  segments. 


d.  Semicircumferential  fliqhts  in  which  point  "B"  is 
throuqh  the  center  of  the  earth  will  also  abort.  Break  such 
segments. 


exactly  opposite  "A" 
flights  into  smaller 


e.  Routes  or  seqments  of  routes  flown  directly  north  or  directly  south  alonq 
a  longitude  line  will  also  abort. 


f.  At  present,  the  model  allows  flights  only  at  altitudes  25,000  feet  and 
35,000  feet. 


g.  All  departure  and  destination  points  must  lie  within  the  latitudinal 
ranqe  from  75°N  to  60°S. 

Restrictions  f  and  a  were  imposed  in  order  to  reduce  the  core  storaae  require¬ 
ments  of  the  model. 
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Chapter  5 


USE  OF  THE  WIND  FACTOR  SIMULATION  MODEL 


5.1  Use  as  a  Subroutine 

The  Wind  Factor  Simulation  Model  is  coded  as  a  constellation  of  seven  FORTRAN 
subprograms  as  described  earlier  in  Table  3.  The  key  subprogram  is  subroutine 
ENRWND,  which  invokes  the  other  six  subprograms  as  necessary.  Thus,  the  wind  factor 
model  is  intended  to  function  as  a  small,  fast  module  within  a  larger  simulation 
model,  such  as  MAC's  COLOSSUS,  as  shown  in  Figure  8. 

The  WFSM  communicates  with  the  user's  simulation  model  only  throuah  the  argument 
list  of  the  user's  call  to  subroutine  ENRWND.*  In  the  user's  model,  wherever  a  wind 
factor  is  desired,  a  call  should  be  placed  to  ENRWND,  supplying  values  for  the  input 
arguments.  Subroutine  ENRWND  then  takes  control,  calculates  a  route-mean  wind 
factor,  and  uses  it  to  prepare  the  output  argument. 

The  user's  call  to  subroutine  ENRWND  should  be  of  the  form, 

CALL  ENRWND  (FRMLAT,  FRMLNG ,  TOLAT,  TOLNG,  JULDAT,  GMT,  PCHRS,  IALT,  IOPTN, 
ASPEED,  GSPEED) 

The  user  selects  calm,  mean,  or  90-percent  worst  winds  by  specifying  a  wind 
option  (IOPTN)  from  0  to  3  (see  Table  6),  and  specifies  the  season  of  the  year  for 
which  he  wants  winds  by  stipulating  a  Julian  date  (JULDAT)  from  1  to  366.  The  user 
then  specifies  the  altitude  of  the  aircraft  by  setting  the  flap  IALT  to  1  for  25,000 
feet  or  to  2  for  35,000  feet.  The  input  variable  ASPEED  indicates  the  reauested 
airspeed  in  knots  and  the  output  variable  GSPEED  indicates  the  calculated  ground 
soeed  in  knots. 


5.2  Inputs 

Input  arguments  for  subroutine  ENRWND  are  presented  in  Table  6.  Special 
formatting  rules  are  as  follows: 

a.  Latitudes  and  longitudes  will  be  in  decimal  degrees,  not  the  degrees- 
minutes-seconds  system. 

b.  South  latitudes  and  east  longitudes  must  be  negative. 


5.3  Output 


Output  arqument  for  subroutine  ENRWND  is  the  variable  GSPEED,  a  real  auantity 
representing  the  aircraft's  ground  speed  in  knots. 


*The  exception  is  that  users  may  call  subroutine  DIfiTAN  directly  if  they  need  a  oreat 
circle  of  distance,  or  function  HDG  if  they  want  an  initial  great  circle  distance. 
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Figure  8.  Hierarchy  of  Models. 

Table  6.  INPUT  ARGUMENTS  TO  SUBROUTINE 

Data 

1 

Name 

Type 

Definition 

Units/Values 

FRMLAT 

Real 

Latitude  of  point  "A" 
(Takeoff) 

Decimal  degrees 

FRMLNG 

Real 

Longitude  of  point  "A" 
(Takeoff) 

Decimal  deqrees 

TOLAT 

Real 

Latitude  of  point  "B" 
(Landing) 

Decimal  deqrees 

TOLNG 

Real 

Longitude  of  point  "B" 
(Landing) 

Decimal  degrees 

JULDAT 

Integer 

Julian  base  date  of 
wind  factor  reouest 

1  -  366 

GMT 

Real 

Greenwich  base  time  of 
wind  factor  reauest 

0.0  -  24.0  hrs 

FCHRS 

Real 

Forecast  hours  ahead 

>0.0  hrs 

IALT 

Integer 

Altitude  index 

1  *  25,000  ft 

2  *  35,000  ft 

IOPTN 

Integer 

Wind  option  or  mode 

0  ■  Calm  wind 

1  »  Simulated  winds 

(defaults  to  3) 

2  »  90-percent  worst 

wind 

3  *  Mean  wind 

AS PE ED 

Real 

Airspeed 

Knots 
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ChaDter  6 

TEST  AND  EVALUATION  OP  THE  WIND  FACTOR  SIMULATION  MODEL 


6 . 1  Original  Tests 

The  Wind  Factor  Simulation  Model  was  tested  over  36  routes.  The  first  22  test 
routes  were  devised  by  the  model  development  team.  They  were  designed  to  test  the 
great  circle  navigation  capability  of  the  model.  Consideration  was  given  to  using 
test  routes  that  would  approximate  the  different  types  of  flights  used  by  COLOSSUS. 
However,  no  attempt  was  made  to  use  actual  airfield  latitudes  or  longitudes.  At 

first,  simple,  within-ouadrant  flights  were  run  from  point  "A"  to  point  "B."  The 

flight  was  then  reversed  and  run  from  "B"  to  "A."  Gradually  the  type  of  flight 

considered  was  made  more  complicated.  The  more  complicated  the  flight  in  terms  of 

distance  and  eouatorial  and  dateline  crossings,  the  more  difficult  it  was  to  solve 
the  great  circle  equation.  Table  7  provides  a  listing  of  all  test  point  coordinates, 
while  Table  8  describes  the  tests  themselves.  Test  routes  are  drawn  in  Figure  9. 

TABLE  7.  TEST  POINT  CORRDINATES, 


LOCATION 

OF 

LOCATION 

OF 

ROUTE 

TEST 

POINT  "A 

II 

POINT  "B 

M 

LENGTH 

NUMBER 

LATITUDE /LONGITUDE 

LATITUDE/LONGITUDE 

(NM) 

1. 

60.  ON 

60. 0W 

30. ON 

100. OW 

2408.7 

2. 

30. ON 

10. OW 

50. ON 

40.  OW 

1807.4 

3. 

50. ON 

10. 0W 

30. ON 

10.  OW 

1807.4 

4. 

10. ON 

20. 0E 

30. ON 

50. OE 

2063.4 

6. 

30.  ON 

50. 0E 

10. ON 

20. OE 

2063.4 

6. 

10. ON 

40. OE 

40. ON 

5.  OE 

2590.6 

7. 

40.  ON 

S.0E 

10. ON 

40.  OE 

2590.6 

R. 

10. OS 

10. OH 

40. OS 

60. OW 

3204.7 

9. 

40. OS 

60.  OW 

10. OS 

10. OW 

3204.7 

10. 

10. OS 

30. 0W 

40. OS 

10. OW 

2093.1 

11. 

40. OS 

10. OW 

10. OS 

30.  OW 

2093.1 

12. 

20. OS 

20. 0W 

50. OS 

70. OE 

2°68. 3 

13. 

50. OS 

70. 0E 

20.  OS 

20. OE 

2968.3 

14. 

50. OS 

30. 0E 

10. OS 

50.  OE 

2599.2 

IS. 

10. OS 

50. 0E 

50. OS 

30.  OE 

2599.2 

16. 

35. ON 

30. 0E 

50.  ON 

10. OE 

1955.5 

17. 

20. OS 

40. 0E 

10. ON 

20.  OW 

3975.5 

18. 

20. OS 

30. 0W 

10. OS 

10. OE 

2449.0 

19. 

20.  ON 

10. OE 

10. OS 

20.  OW 

2627.3 

20. 

10. ON 

170. 0W 

15. OS 

165. OE 

2113.2 

21. 

30. ON 

135. 0E 

45. OS 

120. OW 

7252. R 

22. 

45. OS 

120. 0W 

30. ON 

135. OE 

7252.8 

23. 

46.94N 

67 . 88W( LIZ) 

62. 37N 

0 . 47E { E  GUM ) 

2575.9 

24. 

52.37N 

0. 47E ( EGUN  ) 

46.95N 

67. 88W (LIZ) 

2575.9 

25. 

35.75N 

139.  35E (RJTY  ) 

38.27N 

121 . 93 W ( SUU ) 

4480.0 

26. 

38.27N 

121. 93W{ SSU ) 

35.75N 

139 . 35E ( RJTY ) 

4480.0 

27. 

35.75N 

139.  35E ( RJTY  ) 

13.5RN 

144.92EIPGUA) 

1364.6 

28. 

13.58N 

144. 92E ( PGUA ) 

37.75N 

139.35E (RJTY) 

1364.6 

29. 

41.12N 

?5.90W(0FF) 

35.05N 

106. 62W( ABG) 

623.3 

30. 

35.05N 

106.62W(ABQ) 

41.12N 

95 . 90W ( OFF ) 

62  3.3 

31. 

32.50N 

93 . 67W( BAD) 

41.12N 

1 1 1 . 97W( HI F ) 

1017.5 

32. 

41.12N 

111.97W(HIF) 

32 .50N 

93 ,67W( BAD) 

1017.5 

33. 

35.75N 

1 39 . 35E ( RJTY  ) 

23. 80S 

133. 90E ( ASI S ) 

3588.7 

34. 

23. 80S 

1 33 . 90E { ASI S  ) 

35.75N 

139 . 35E ( RJTY ) 

3588.  7 

35. 

63.97N 

22 . 60W  ( B I K F ) 

61.30N 

149 . 80W ( EOF ) 

2920.5 

36. 

61.30N 

149.80W(E0F) 

63.97N 

22 . 60W ( B I K F  ) 

2920.5 

TABLE  8.  WIND  FACTOR  SIMULATION  MODEL  TESTS 


TESTING 

SESSION  TEST 

NUM 

1  2  3  4  5  6  7  R  9  10  11  12 

FSSSSSSFFSS 
FSSSSSSSSSS 
SSSSSSSSSSS 
SSSSSSSSSSS 
SSSSSSSSSSS 
SSSSSSSSSSS 

13  14  15  16  17  18  19  20  21  22  23 

A  SSSSFSSSFFS 

B  SSSSSSSSFFS 

C  SSSSSSSSFFS 

D  SSSSSSSSFFS 

E  SSSSSSSSSSS 

F  SSSSSSSSSSS 

25  26  27  28  29  30  31  32  33  34  35 


A 

R 

C 

0 

E 

F 


S  *  SUCCESSFUL  TEST 
F  *  FAILED  TEST 
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6.2  Additional  Tests 

It  was  the  great  circle  computations  (specifically  the  first  guess  longitude  and 
the  longitude  equator  crossing  calculations)  that  slowed  progress  on  completing  the 
model.  In  addition,  14  more  test  runs  were  recommended  by  the  Aerospace  Sciences 
Branch  of  USAFETAC.  These  new  tests  were  different  from  the  original  22  in  that 
soecific  COLOSSUS  terminals  were  used  as  points  "A"  and  "B."  While  the  oriqinal  22 
tests  were  successful  (20  of  22  were  successful  the  first  week;  see  Table  7  and 
Figure  9),  the  new  tests  posed  a  greater  challenge.  In  fact,  seven  of  the  14  failed 
the  first  time  they  were  run. 

These  failures  caused  problems  not  iust  with  the  specific  run  that  failed,  but 
also  with  earlier  successful  runs.  Debugging  a  particular  test  run  often 
inadvertently  caused  previously  successful  runs  to  fail.  Consequently,  the 
development  team  had  to  test  every  one  of  the  36  runs  whenever  any  changes  were  made 
in  the  wind  factor  program  code. 

Midway  throuqh  the  test  procedure,  the  wind  factor  module  passed  32  of  the  36 
test  runs.  Two  of  the  unsuccessful  runs  were  test  21,  (30.0°N,  135. 0°E)  to  (45.0°S, 
120. 0°w),  and  its  reverse,  test  22.  The  other  two  unsuccessful  tests  were  numbers  33 
and  34.  Tests  33  and  34  are  the  Yokota  AB,  Japan  (35.8°N,  139. 4°E)  to  Alice  Sprinqs, 
Australia  (23.8°S,  133. 9°E)  route  and  its  reverse.  Test  21/22  was  a  very  long  route, 
over  7000  miles.  The  route  was  further  complicated  by  it's  crossing  both  the  equator 
and  the  dateline.  Test  33/34  was  complicated  in  that  it,  too,  was  lonq  and  crossed 
the  equator,  but  more  so  in  that  it  was  almost  a  due  north-south  route.  In  fact,  the 
3588  miles  were  traversed  with  a  lonaitudinal  variation  of  less  than  6  degrees.  This 
lack  of  longitudinal  variation  approaches  one  of  the  limitations  in  the  wind  factor 
module,  i.e.,  a  direct  north-south  route  is  not  allowed  and  will  abort  the  proaram. 

6 . 3  Problems  Encountered 


In  these  more  complicated  tests,  flights  proceeded  normally  until  they  reached 
the  equator.  At  that  point  the  aircraft  appeared  to  "jump"  from  the  equator  to  point 
"B“  ignorinq  the  intervening  or  id  crossings.  For  example,  the  flight  from  Yokota  to 
Alice  Springs  (test  33)  picked  up  the  points  from  Yokota  southward  to  the  equator  but 
then  "jumped"  from  the  equator  to  Alice  Springs.  Thus,  it  picked  ud  points  at 
Yokota,  (30.0°N,  138. 7°E),  (15.0°N,  137. 2°E)  and  the  equator  crossing  (0.0°, 
136. 0°F).  It  then  "jumped"  from  the  equator  crossing  to  Alice  Sorings,  missing  the 
point  at  (15.0°S,  134. 7°E).  When  the  route  was  reversed  (test  34),  the  module  picked 
up  the  points  from  Alice  Springs  northward  to  (15.0°S,  13A.7°E)  and  then  to  the 
equator  crossing  where  it  "jumped"  to  Yokota,  thus  missing  the  points  at  (15.0°N, 
137. 2°E)  and  (30.0°N,  138. 7°E).  Since  these  points  are  actually  grid  crossinqs  or 
crossings  of  boundaries  of  the  wind  factor  grid  sectors,  their  mean  wind 
contributions  to  the  calculation  of  the  overall  wind  factor  were  beina  ignored.  Even 
though  a  route  wind  factor  was  produced  it  was  in  error. 

These  errors  were  corrected  by  reversing  the  direction  of  the  spherical  triangle 
as  the  aircraft  crossed  the  eouator.  This  change  solved  the  problem,  and  all  four  of 
these  test  runs  were  successful.  However,  when  the  entire  36  tests  were  run,  test  35 
(Keflavik,  Iceland,  to  Elmendorf,  Alaska)  failed.  The  problem  here  involved  a  bad 
computation  of  the  fictitious  equator  crossing.  Once  the  problem  was  isolated  and 
corrected,  all  36  tests  were  run  again,  this  time  successfully. 
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Chapter  7 
APPLICATIONS 


7.1  Use  in  COLOSSUS 

The  Wind  Factor  Simulation  Model  was  designed  to  a  response  to  a  simulation 
requirement  of  the  Military  Airlift  Command  (MAC).  The  wind  factor  module  was  needed 
as  an  element  of  the  Military  Airlaft  Command  Resource  Optimization  (MACRO)  M-14 
simulation  project  called  COLOSSUS,  which  attempts  to  simulate  the  capability  of 
Military  Airlift  Forces  to  respond  to  a  real  world  contingency  anywhere,  anytime. 

7.2  Potential  Use  in  Other  Simulators 


Although  designed  in  response  to  MAC'S  COLOSSUS  requirement,  the  WFSM  is  oeneral 
enough  to  make  it  useful  for  other  purposes  and  flexible  enouqh  to  give  it  growth 
potential.  Other  simulators  that  require  wind  factor  information  mav  be  able  to  make 
use  of  the  present  model  quite  easily.  As  resources  for  conducting  real  world  tests 
and  exercises  gradually  decrease,  the  advantages  of  simulation  qrow  more  apparent. 
Since  the  very  concept  of  simulation  is  to  model  real  world  events  representatively, 
it  follows  that  including  weather  in  simulation  of  environmentally  sensitive  Air 
Force  operations  is  advisable.  As  Air  Force  simulation  models  become  more  complex 
and  more  representative  of  the  real  world,  requirements  such  as  that  of  COLOSSUS  will 
increase  in  number  and  grow  in  complexity. 

7.3  In-house  Fast  Response  wind  Factor  Capability 

USAFETAC  receives  many  requests  for  flight  planning  information  and  wind  factors. 
At  present  an  in-house  computer  program  is  used  to  respond  to  these  requests,  while 
that  method  is  proven  and  has  been  in  use  for  a  couple  of  years,  it  is  costly  in 
terms  of  computer  run  time  and  core  storage.  With  improvements  in  resolution,  the 
WFSM  may  prove  to  be  an  economical  replacement  in  many  cases  for  the  existing  wind 
factor  program. 
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Appendix  A 

SPHERICAL  AND  GLOBAL  COORDINATE  SYSTEMS 


A  special  orthogonal  curvilinear  coordinate  system,  the  spherical  coordinate 
system  uses  coordinates 
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The  spherical  coordinates  (r,  0  ,  cp)  are  related  to  global  coordinates  (a  +  z, 
LAT,  LNG) ,  where 
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To  convert  from  global  coordinates  to  spherical  coordinates,  eauation  set  (A-l) 
is  used.  To  convert  from  spherical  coordinates  to  global  coordinates,  eauation  set 
(A-2)  is  used.  Since  r  ■«  a,  the  conversions  are  between  0  and  latitude,  and  cp  and 
longitude. 

Spherical  coordinate  cp  is  the  azimuth  angle  and  is  measured  in  radians  counter¬ 
clockwise  from  Greenwich  looking  down  from  space  onto  the  North  Pole.  The  absolute 
value  ranges  from  0  to  2  tt  radians.  Negative  values  Indicate  clockwise.  (See  Picture 
A-l) . 
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Appendix  A 
THE  BETA  ANGLE 


The  meteorological  wind  direction  dd  can  be  expressed  in  the  form  of  a  beta 
anqle,  g  .  The  beta  anqle  describes  the  direction  toward  which  the  wind  is  blowing 
(heading)  in  terms  of  radians  counterclockwise  from  eastward  (96°)  . 

This  system  is  different  from  the  conventional  method  of  indicating  wind 
direction.  The  usual  manner  is  to  express  wind  as  a  direction  from  which  the  wind  is 
blowing.  The  advantage  in  using  the  beta  angle  scheme  is  that  calculations  involving 
aircraft  heading  and  beta  angle  wind  direction  (a  heading  direction)  are  more  easily 
handled  when  the  frame  of  reference  is  directionally  the  same. 

The  following  illustration  serves  to  show  the  relationship  between  conventional 
wind  direction  and  the  beta  angle  method  of  indicating  wind  direction. 

Direction  Beta  Angle 

Toward  Which  (degrees)  Beta  Angle 

(0°  =  N,  90°  =  E,  etc.)  (CCW  from  Eastward)  (Radians) 


LIST  OF  ABBREVIATIONS  AND  ACRONYMS 


AS PEED 

BRLNG 

BRLAT 

COLOSSUS 

COMMON 

CPU 

01 STAN 

DIRM 

ENRWND 

FCHRS 

FORTRAN 

FRMLAT 

FRMLNG 

GCD 

GRTCIR 

GSPEFD 

GMT 

HDG 

HDGM 

HP 

IALT 

ISEASN 

IOPTN 

JULDAT 

MAC 

MACRO 

NM 

NRND 

NSEC 

STOP 

SPHGLO 


airspeed 

FORTRAN  suborogram 
FORTRAN  subprogram 
MAC  airlift  system  simulation 

Area  of  computer  storage  shared  between  the  main  program  and  its 
subprograms 

Central  processing  unit 
FORTRAN  subprogram 
wind  direction 
FORTRAN  subprogram 
forecast  hours  ahead 

Formula  Translation  Language,  a  computer  language  suitable  for 
scientific  problem  solving 

from  latitude 

from  longitude 

great  circle  distance 

FORTRAN  subprogram 

ground  speed 

Greenwich  Mean  Time 

FORTRAN  subprogram 

aircraft  heading  at  location  (XLATM,  XLNGM) 

Hewlett-Packard  Company 
Altitude  index 
season  index 
wind  option 
Julian  Date 

Military  Airlift  Command 

Military  Airlift  Command  Resources  Optimization 

Nautical  miles 

record  number 

grid  sector  number 

Single  Integrated  Operation  Plan 

FORTRAN  subprogram 
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SNP 

SPDM 

TOLAT 

TOLNG 

USAFETAC 

VARM 

VGM 

VCM 

VGACC 

VCACC 

VACC 

WFSM 

WBAR 

WBARBR 

XLATC 

XLNGC 

XLATNP 

XLNGAS 

XLATNX 

XLNGNX 

XLATAS 

XLATM 

XLNGM 


distance  between  Current-Point  and  Next-Point 
wind  speed 
to  latitude 
to  longitude 

United  States  Air  Force  Environmental  Technical  Applications  Center 
wind  variance  at  location  (XLATM,  XLNGM) 

along-oround-track  wind  component  at  location  (XLATM,  XLNGM) 

cross-ground-track  wind  component  at  location  (XLATM,  XLNGM) 

sum  of  along-track  component  for  leg  times  leg  length 

sum  of  cross-track  component  for  leq  times  leg  length 

sum  of  wind  variance  times  leg  lenqth 

Wind  Factor  Simulation  Model 

time-averaged  wind  factor 

route-mean  wind  factor 

Latitude  of  Current-Point 

Lonqitude  of  Current-Point 

Latitude  of  Next-Point 

Longitude  associated  with  XLATNX 

Latitude  of  Candidate  #1  for  Next-Point  (next  latitude  crossina) 
Longitude  of  Candidate  12  for  Next-Point  (next  longitude  crossing) 
Latitude  associated  with  XLNGNX 
Latitude  midpoint 
Lonqitude  midpoint 
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